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Abstract 

We investigate quantum tunneling through a single electric and/or magnetic barrier on the surface of a three- 
dimensional topological insulator. We found that (1) the propagating behavior of electrons in such system exhibits 
a strong dependence on the direction of the incident electron wavevector and incident energy, giving the 
possibility to construct a wave vector and/or energy filter; (2) the spin orientation can be tuned by changing the 
magnetic barrier structure as well as the incident angles and energies. 
PACS numbers: 72.25.Dc; 73.20.-r; 73.23.-b; 75.70.-i. 



1. Introduction 

The recent discovery of a new quantum state of matter, 
topological insulator, has generated a lot of interest due 
to its great scientific and technological importance [1-5]. 
In a topological insulator, spin-orbit coupling opens an 
energy gap in the bulk, and results in helical surface 
states residing in the bulk gap in the absence of mag- 
netic fields. Such surface states are spin-dependent and 
are topologically protected by time-reversal symmetry 
[4-7] and distinct from conventional surface states, 
which are fragile and depend sensitively on the details 
of the surface geometry and bonding. This discovery 
sparked intensive experimental and theoretical interests, 
both for its fundamental novel electronics properties as 
well as possible applications in a new generation of elec- 
tric devices. 

Very recently, the surface states in Bi-based alloys, Bi^ 
^Sb^, Bi 2 Se 3 , Bi 2 Te 3 , were theoretically predicted [6,8] 
and experimentally observed by using angle-resolved 
photoe-mission spectroscopy (ARPES) [9-12]. These 3D 
topological insulators have robust and simple surface 
states consisting of a single Dirac cone at the T point 
[8]. Note that this Hamiltonian appears similar to gra- 
phene [13], but topological insulators have an odd num- 
ber of massless Dirac cones on the surface, ensured by 
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the Z 2 topological invariant of the bulk, while graphene 
has twofold massless Dirac cones at the K and IC valleys. 
Another essential difference is connected with spin- 
related properties. In the surface, Hamiltonian of the 3D 
TI a acts on the real spin of the charge carriers, while 
for graphene it stands for the pseudo spin, i.e., the A 
and B sublattices of graphene. Hence, it is natural to 
manipulate spin transport on the surface of a 3D topo- 
logical insulator by controlling the electron orbital 
motion. Based on the topological surface Hamiltonian, it 
is clear that a • k is a quantum conserved quantity 
which implies that spin and momentum of the electron 
are locked. For instance in a tunneling process, the 
reflected electron will reverse its spin due to the helical 
property of the surface states, i.e., the spin-momentum 
locking [14]. This feature will lead to some interesting 
phenomena, such as the spin-dependent conductance 
[15] and the twisted RKKY interaction [16]. 

In this study, we investigate electron tunneling 
through single electric and magnetic potential barriers 
which can be created by depositing a ferromagnetic 
metallic strip on the surface of a 3D topological insula- 
tor. We find that the in-plane spin orientation of the 
transmitted and the reflected electrons can be rotated 
over certain angles that are determined by the incident 
angle and energy. Our results demonstrate that the mag- 
netic field of the magnetic barrier bends the trajectory 
of the electrons, and therefore rotate the spin. 
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2. Theory 

First we focus on the electron transmission through a 
single electric and magnetic potential barrier on the sur- 
face of a 3D topological insulator, as shown in Figure 1. 
We can create magnetic and electric potentials under- 
neath a superconducting plate and a ferromagnetic strip 
[17-20]. The relevant magnetic and electric fields are 
directed perpendicular to the surface of the 3D TI, i.e., 
V(x) = V[®(x) - ®(x - d)]/2 and B z (x) = B[®(x) - ®(x - 
d)\\% where 0 is the Heaviside step function. The mag- 
netic field is perpendicular to the surface and the 
Landau gauge A = (0, Bx, 0) is adopted in our calcula- 
tion. The low-energy electrons near the T point of the 
Dirac cones can be well described by the effective 
Hamiltonian [21], 



H = v?o • (n xz) + V + H z 



(i) 



where v F is the Fermi velocity, &{i = x, y, z) are the 
Pauli matrices, V is the gate voltage applied on the mag- 
netic metal strips, and the last term H z = gfi B <j - B is 
induced by Zeeman spin spitting. Note that for g = 23 
in Bi 2 Se 3 , the Zeeman term affects the transmission 
slightly at low magnetic field. The momentum is 7i = p 
+ eA, where the vector potential of the inhomogeneous 
magnetic field generated by the magnetic metal stripe, 
z^, is the unit vector normal to the surface. Note that 
[a • (n x z),H] = 0, which implies that a • (n x z) is a 
quantum conserved quantity during the tunneling pro- 
cesses, i.e., spin-momentum locking. If the incident elec- 
trons are spin polarized along the direction of the 
vector (tt x z), the magnetic field bends the trajectory of 
the electrons, resulting in a rotation of the spin of the 
transmitted electrons. The reflected electrons suffer 
similar spin rotations accompanied by the reversal of 
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Figure 1 Schematic of an electric and magnetic potential 
barrier on the surface of a 3D TI. 



the momentum p x . Interestingly, the gate voltage can be 
used to control the reflection and transmission, and 
therefore tune the spin polarization of the reflected and 
transmitted electrons. 

For simplicity, we introduce the dimensionless units: 
l B = [hleB 0 ] 112 , B(x) -> BqB(x), A y (x) -> BqIbA^x), E -> 
EqE, r — » l B r, k — » k/l B , and rewrite the Hamiltonian as 



H : 



V + gfl B B ky+Ay 



ky+Ay 



ik x V 



+ i% 
g^BB 



) 



(2) 



In the presence of the magnetic-electric potentials, the 
wavevector satisfies 
k 2 x + (ky + A y ) 2 = (E — V) 2 - (g/x B B) 2 . Note that the 
translational invariance along the y direction gives rise 
to conservation of k y) and thus the solutions can be 
written as ^/(x,y) = x//[x)e ik y y ^ In the free region, the 
eigenvalues are E ± = ±v F n, and the corresponding eigen- 
vectors are 



f L {x) = 



1 
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In the barrier region the vector potential is A = (0, Bx, 
0) and the solution can be expressed in terms of para- 
bolic cylinder functions D v , 



D v /2-i[±y/2(x + k y )] 
±J±D v/2 [±V2(x + k y )] 



(5) 



with complex coefficients c+ and v = (E F - V) 2 - (gfi B B) 

2 . After some lengthy algebra, all of the above coeffi- 
cients of the wave functions can be obtained from the 
boundary conditions. 

3. Spin and momentum filtering 

First we consider that the incident electrons are spin 
polarized along the direction of the vector (it x z), i.e., 
perpendicular to the direction of the electron motion. 
We examine the transmission characteristics for a pure 
electric barrier (B = 0, V = 4), which are shown in Fig- 
ure 2a,b. In Figure 2a we plot the transmission probabil- 
ity as a function of the incident angle and energy. The 
symmetry of the low-energy effective Hamiltonian when 
A = 0 ensures the invariance of the transmission with 
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Figure 2 Transmission spectrum through the single magnetic/ 
electric barrier and spin rotation of transmitted/reflected 
electrons, (a) The contour plot of the transmission probability as a 
function of the incident angle and the incident energy, for a fixed 
barrier width d = 1 and height V=4, for magnetic field 5 = 0. (b) The 
spin orientation of the incident (black arrows), transmitted (red 
arrows), and reflected (green arrows) electrons as a function of the 
incident angle and energy. The lengths of the arrows indicate the 
transmitted or reflected probabilities. The magnetic field unit is B 0 = 
1 T, the energy unit is E 0 = 16 meV, and the length unit is l B = 26 nm. 



respect to k y — > -k y . Note that, perfect transmission 
always exists in the vicinity of normal incidence, i.e., 
this is the so-called Klein tunneling, induced by the heli- 
cal property of the Dirac fermion. This perfect tunneling 
process can even occur for a low incident energy, and a 
high and wide barrier. Figure 2b shows the spin orienta- 
tion as a function of the incident angles and the incident 
energies. From the spin orientation one can see that the 
transmitted electron spins are polarized along the same 
direction as the incident electron spin, indicating that a 
pure electric barrier without magnetic field will not 
affect the spin orientation, the reason is that an electric 
barrier cannot bend the trajectories of the electrons and 
the outgoing electrons will propagate in the same direc- 
tion as the incident ones. The spin of the reflected elec- 
tron, with reversed longitudinal wavevector and 
conserved transverse wavevector, is rotated due to the 
spin-momentum locking. 



Next we consider the tunneling process through a pure 
magnetic single barrier. The transmission probability is 
shown in Figure 3a for a magnetic barrier (B = 1) with 
width d = 1. Compared to the pure electric barrier case 
(see Figure 2a), the transmission becomes asymmetric with 
respect to the in-plane momentum k y parallel along the 
interface, since the magnetic field breaks the time reversal 
symmetry. Note that Figure 3a implies that for a certain 
incident angle (p, tunneling is forbidden, the boundary of 
the total reflection region (T = 0) can be approximately 
given by the relation k y + Bd = E, which implies that the 
transmitted wave is an evanescent mode which decays 
exponentially along the propagating direction. This total 
reflection is always present regardless of the incidence 
angle (p when E <Bd\% and results in an asymmetric beha- 
vior of the transmission as a function of the incident angle. 
This total reflection can be understood semiclassically: 
when the width d of the barrier is large as compared to the 
Fermi wavelength X F , the electron inside the barrier will 
move on a cyclotron orbit, and if the cyclotron orbit radius 
R c <d, the incident electron will exit the barrier region 
backwards eventually. This feature illustrates that the Dirac 
quasi-particles can be confined by the inhomogeneous 
magnetic field, and, more interestingly, the incident elec- 
tron spin-oriented parallel along the interface will be 
flipped in the total reflection region. Figure 3b shows the 
spin orientation as a function of the incident angles and 
energies, which is also asymmetric with respect to the inci- 
dent angles, as a consequence of the perpendicular inho- 
mogeneous magnetic field. Note that the transmitted 
electron spin no longer points along the same direction of 
the incident electron spin, but is rotated over an angle, 
which is determined by the incident angle, the electron 
energy and the barrier width. In addition, changing the 
length of the barrier and/or the magnetic field can tune the 
total reflection and the perfect transmission regions. 

We also examined the electron transmission through a 
combined electric and magnetic barrier. The transmission 
and spin orientation (see Figure 4a,b) become very differ- 
ent from that of the pure magnetic barrier case. The inter- 
play between electric barrier and magnetic field strongly 
reduces the perfect transmission region, and provides us 
with an additional way to control the transmission and the 
spin orientation of the transmitted and reflected electrons. 

The transmission feature shown above can be exam- 
ined by the measurable quantities, the conductance G, 
and Fano factor F [22,23]. The ballistic conductance and 
Fano factor for a given Fermi energy at zero tempera- 
ture are given by summing over the modes, 



G(E F ) = GoV°° T„(E F ), 

- J n=— oo 

F(B F ) = T »( £ f)[1 " T„(E F )]/ T„(£ F ), 

-^n=— oo *—^n=—oo 



(6) 
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Figure 3 The same as Figure 2, but now for a pure magnetic 
barrier with 8 = 1, and V = 0. 
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Figure 4 The same as Figure 2, but now for a combined 
electric and magnetic barrier with 8 = 1, and V = 4. 
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Figure 5 Conductance and Fano factor of the single magnetic/electric barrier, (a) The conductance as a function of the incident energy 
with B = 1, V = 0 (black solid line); B = 1, V = 4 (red solid line); and B = 0, V = 4 (blue solid line), (b) The Fano factor as a function of the 
incident energy for the same magnetic/electric barrier. 
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where G 0 = — is taken as the conductance unit. The 
7th 

dependence of the conductance G and the Fano factor F 
on the incident energy E ¥ are plotted in Figure 5. In 
low-energy region, the magnetic barrier can suppress 
the electron transmission probability for any incident 
angles, and thus lead to vanished conductances. The 
conductance increases as the Fermi energy since the 
magnetic field is not strong enough to confine the elec- 
trons. For a pure electric barrier, the conductance is 
very large in low energy region, since the electrons can 
effectively tunnel through the barrier via the down 
branch of the surface state, i.e., the Klein tunneling. 
When the incident energy approaches the top of an 
electric barrier, the conductance decrease drastically. 
This is because the wave vectors in the electric barrier 
become imaginary if the electrons are not normal inci- 
dent. The Fano factor shows good correspondence with 
the changes of conductance. We can see clearly a signif- 
icant change when the barrier changes from transparent 
to opaque. These quantities are what experimentalists 
can actually measure. 

4. Conclusion 

In summary, we investigated theoretically quantum tun- 
neling processes through a single electric and magnetic 
barrier on the surface of a 3D topological insulator. Our 
theoretical results show that the propagating behavior of 
electrons in such a 2D system can be controlled by elec- 
tric and magnetic barriers. We have also observed the 
rotation of the electron spin, depending on the para- 
meters of the magnetic barrier structure as well as the 
incident angle and energy. This investigation could be 
helpful to offer a promising functional unit for future 
wavevector, and/or spin filtering quantum devices. 
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